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The principle part of Einstein equations in the harmonic gauge consists of a constrained system 
of 10 curved space wave equations for the components of the space-time metric. A well-posed initial 
boundary value problem based upon a new formulation of constraint-preserving boundary conditions 
of the Sommerfeld type has recently been established for such systems. In this paper these boundary 
conditions are recast in a geometric form. This serves as a first step toward their application to 
other metric formulations of Einstein's equations. 
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It is extremely satisfying to contribute this article in acknowledgment of the important influence that Jiirgen Ehlers 
has had on my career. When I was a first year graduate student at Syracuse University in the Fall of 1959, my PhD 
adviser Peter Bergmann had collected an astonishing percentage of the young talent engaged in general relativity. 
Besides the tenured faculty members Art Komar and Richard Arnowitt, the postdocs and visiting scientists included 
Roy Kerr, Wolfgang Kundt, Ted Newman, Roger Penrose, Asher Peres, Ivor Robinson, Englebert Schucking, Andrej 
Trautman and Manfried Trumper, as well as Jiirgen. At the end of that first year, when Bergmann called me to his 
office to discuss thesis research topics, from this group of experts it was Jiirgen whom he invited to join in offering me 
advice. In response to Bergmann's opening question "So, what would you like to work on?" , I suggested singularity 
structure. I had recently enjoyed reading the nice paper by Bergmann and his former student Ray Sachs on singularities 
in linearized gravitational theory, as well as the Einstein-Infeld- Hoffman paper on the motion of singularities. Jiirgen 
immediately took over the conversation, "What is your underlying motivation?" I naively offered the possibility that 
singularity structure might be used to understand elementary particles. I then received a wide ranging lecture from 
Jiirgen who carefully explained the conceptual problems underlying singularities even in electromagnetic theory and 
why the situation was much more complicated in general relativity. At that time, a global picture of the structure of 
the Schwarzschild singularity, which could not be anticipated from linearized theory, had just emerged with the help 
of Kruskal's extension of the exterior spacetime. The geometric approach of classifying singularities in terms of the 
incompleteness of the spacetime manifold had just begun. The important distinction between spacelike, timelike or 
null singularities was not yet recognized. In the following years, an effective approach to this difficult subject slowly 
developed from a great deal of effort by many people. Fortunately for me, I accepted Jiirgen's impromptu lecture as 
good advice to steer clear of the subject in my graduate research. 

Also at that time, general relativity was a small field going through a renaissance, which was centered in the 
United States about Bergmann at Syracuse and John Wheeler at Princeton. Bergmann's research, which was focused 
on the quantization of gravity, had led to several reformulations of Einstein's theory in terms of Lagrangian, phase 
space and Hamilton- Jacobi methods. But in all of these formalisms the same problems associated with dealing with 
nonlinearity, identifying the proper observables and handling the constraints posed the same stumbling blocks against 
a real physical understanding of quantum gravity that persist today. However, there were several other current 
developments which would have seminal impact on the future of the field: the Kruskal extension of Schwarzschild 
spacetime, the Kerr solution, the characteristic description of gravitational radiation by Bondi and Sachs, followed by 
Penrose's conformal version; Penrose's spinor description of gravity, which later led to twistor theory; and Jiirgen's 
geometric reformulation of general relativistic hydrodynamics and thermodynamics, which set the standard for the 
ensuing transition in general relativity from the ambiguities of a coordinate dependent, calculational approach to 
a geometric approach. The discovery of quasars in the following years would bring these developments in general 
relativity to the attention of astrophysicists, who would become a captive audience to lectures on the geometry of 
curved spacetimes, thanks to Jiirgen's help in organizing the first "Texas Meeting" . 

Throughout my career I have benefited much from Jiirgen's advice, especially during my two sabbaticals with 
his group in Garching and my annual visits to the Albert Einstein Institute, an institute which was created due to 
his efforts as founding director. Through my exposure to his standards of clarity and rigor and his emphasis on a 
geometric picture, Jiirgen endures as my mentor. I know that he would much prefer to be the cause of an article 
of substance rather than praise, so in that spirit I present the following discussion of the geometrical aspects of the 
initial-boundary value problem (IBVP) in general relativity. Quite some time ago, Jiirgen and his student Saskia 
Kind [l[ treated the spherically symmetric IBVP for a general relativistic fluid, and they later applied the work to 
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stellar oscillations in collaboration with Bernd Schmidt [2|. At that time, little was known about the treatment of 
boundaries in general relativity outside of the spherically symmetric case. Here I treat the vacuum problem in the 
absence of symmetries. In recent work, catalyzed by interactions at the Albert Einstein Institute, the well-posedness 
of the IB VP for Einstein's equations has been established using harmonic coordinates [|| . This puts the IB VP on the 
same analytic footing as the Cauchy problem, whose well-posedness was also established using harmonic coordinates in 
the classic work of Choquet-Bruhat [4|]. However, the geometric formulation of the boundary conditions and boundary 
data for the IBVP is more complicated than for the Cauchy problem, in which the initial data can be expressed in 
terms of the intrinsic metric and extrinsic curvature of the initial Cauchy hypersurface. In this article I will discuss 
these geometrical considerations which underlie the treatment of a boundary. 

In the Cauchy problem, initial data on a spacelike hypersurface So are extended to a solution in the domain of 
dependence T>(So) (which consists of those points whose past directed characteristics all intersect So)- In the IBVP, 
data on timclike boundary T transverse to So is used to further extend the solution to the domain of dependence 
V(S UT). 

The IBVP for Einstein's equations has only relatively recently widespread attention, when its importance to numer- 
ical relativity was pointed out [f| ■ The first well-posed IBVP was achieved for a tetrad formulation of Einstein's theory 
using a first differential order system which included the tetrad, the connection and the curvature tensor as evolution 
fields Q . Subsequently, a well-posed IBVP was formulated for the harmonic formulation of Einstein's equations as 
second order wave equations for the metric 0] . This extended the classic analytic treatment of Choquet-Bruhat Q to 
the well-posedness of the harmonic IBVP. The initial data for the Cauchy problem has a simple description in terms 
of the intrinsic metric and extrinsic curvature of So. The aim of this article is to present an interpretation of the 
boundary data in the IBVP in terms of the geometry of T. This is not only important for a clearer understanding 
of the nature of the gravitational IBVP but also, from a practical point of view, for the application of boundary 
conditions in numerical relativity. The boundary conditions developed in Q have been successfully implemented in 
a evolution code based upon the harmonic formulation of Einstein's equations 0]. However, much of the numerical 
work in general relativity is carried out using other metric formulations, e.g. the BSSN formulation [1, where 
well-posedness of the IBVP currently remains an unresolved issue. The geometric formulation of boundary conditions 
for the metric presented here is a step in that direction. 

I begin with a short review of the Cauchy problem in SecJH followed by a discussion of the new difficulties presented 
by the IBVP and their analytic resolution using a harmonic formulation in See's. [Hi and Iffll where the analytic form 
of the boundary conditions is expressed in terms of partial derivatives of the harmonic metric. In See's. HVl and IVl 
these are recast in covariant form in terms of geometric structures intrinsic to the IBVP. 

I. THE CAUCHY PROBLEM 

The initial data for the Cauchy problem can be formulated in a purely 3-dimensional form in terms of the intrinsic 
metric /i M „ and extrinsic curvature fc^ t „ of the initial Cauchy hypersurface. Here, for notational simplicity, I will use 
Greek indices rather loosely to describe either 3-dimensional or 4-dimensional objects. A major notational complication 
of the IBVP arises from the 3 + 1 decomposition intrinsic to the Cauchy hypersurfaces and the separate 3 + 1 
decomposition intrinsic to the timelike boundary T . As consistently as possible, I will use lower case letters, e.g. h^ v 
and kpu, for geometric objects associated with the Cauchy hypersurfaces and upper case letters for their counterpart 
on the boundary. The Cauchy data are subject to the Hamiltonian and momentum constraints 

= ^R+(k£f~ [=2G^n u ) (1.1) 

and 

= ^(A? - (= (1.2) 

where is the covariant derivative and is the curvature scalar associated with h^ v . Subject to these 

constraints, the Cauchy data determine a solution of Einstein's equations which is unique up to a diffcomorphism 
(cf. [l(| for an exposition with many techniques common to the approach adopted here for the IBVP.) 

However, as already hinted by the parenthetical appearance of the Einstein tensor in (jl.l[) and (|1.2|l . this 
disembodied 3-dimensional form of the Cauchy data hides the complexity of the underlying 4-dimensional space-time 
problem. In order to evolve the data it is necessary to introduce a foliation of the spacetime by Cauchy hypersurfaces 
St, with unit timelike normal n^. The evolution of the spacetime metric is carried out along the flow of a timelike 
vector field f related to the normal by the lapse a and shift (3^ according to 



= an» + , /Fr^ = 0. 
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In numerical applications, the evolution is coordinatized by a time function t satisfying Ctt = 1 and — — aV M f, 
and spatial coordinates satisfying £ t x l — 0, where C t is the Lie derivative with respect to t^. In this 4-dimensional 
setting, h liv n 1 ' = and the spacetime metric is given by 

An additional complexity, already apparent from the arbitrariness of the diffeomorphism freedom in the solution, 
is that Einstein equations do not directly provide a hyperbolic system of evolution equations. For that purpose, it 
is necessary to restrict the gauge freedom. This can be accomplished by introducing harmonic coordinates, i.e. four 
solutions x^ = (t,x l ) of the curved space scalar wave equation DgX^ = 0, so that in these coordinates the metric 
satisfies 

C := = gP°Y% = -_ -Lft,^ = 0, (1.3) 



where 7^ = \]—gg^ v - 

This leads to the standard harmonic reduction of the Einstein tensor 

EilV ._ qhu _ w (n T v) + l g »«V p TP, (1.4) 

whereby the harmonic conditions (| 1 - 3[) together with the vacuum Einstein equations give rise to the quasilinear system 
of coupled wave equations 

= 2j—gE» v = g pa d p d a ^ v + N^, (1.5) 



where represents terms which do not enter the principle part. 

In the harmonic formulation, the Hamiltonian constraints are satisfied as a consequence of the harmonic constraints 
(|1.3[) provided the initial data satisfies r M | t= o = and 8tT^\t=o — 0. The harmonic conditions also determine dta 
and <9t/3, so that the remaining freedom in the initial gauge data reduce to a|t = o and /3 M |t = o We return to the issue 
of constraint preservation in Sec. IVIII 

In summary, the Cauchy data necessary for determination of a unique space-time metric consist of g^, v \t=o — and 
dtg^v\t=o — 0, subject to constraints. This 4-dimensional space-time version of the initial data is referred to as the 
thin sandwich formulation, as opposed to the disembodied 3-dimensional version in which h^lt—o = k^\ t= o — are 
prescribed. The resulting harmonic evolution is not only unique but the solution depends continuously on the choice 
of initial data, i.e. the harmonic Cauchy problem is well-posed jij. 



II. DIFFICULTIES OF THE BOUNDARY TREATMENT 



The IBVP has quite different features than the Cauchy problem, as can be inferred from the simple properties of 
the flat-space scalar wave equation 

{8^-5^8,8^ = x l = (x,y,z) 

in the region x < 0. For a given propagation direction k 1 , there will be two characteristics (light rays) x % = ±kH 
crossing the Cauchy hypersurface at t = 0, but only one characteristic with k x > crossing the boundary at x = 0. 
As a result, although the initial Cauchy data consist of <I> | * = o and dt$\t=o, only half as much boundary data can be 
freely prescribed at x = 0, e.g the Dirichlet data dt$\ x =o, or the Neumann data <9 2 ,<I , | ;r= o or the Sommerfeld data 
(8 t + 8 x )&\ x =o (based upon the derivative in the outgoing characteristic direction). In the gravitational case, this 
inability to prescribe both the metric and its normal derivative on a timelike boundary implies that you cannot freely 
prescribe both the intrinsic metric of the boundary and its extrinsic curvature. In terms of the metric components, 
the most you can describe is a single quantity, e.g. the Sommerfeld data K^d^gpc where if M is an outgoing null 
direction. Such a Sommerfeld boundary condition is most beneficial for numerical work since it allows discretization 
error to propagate across the boundary (whereas Dirichlet and Neumann boundary conditions reflect the error and 
trap it in the grid). 

Inability to prescribe both the metric and its normal derivative, complicates constraint enforcement on the boundary, 
i.e. the Hamiltonian and momentum constraints cannot be enforced directly because they couple the metric and its 
normal derivative. Instead, the approach in [3] is to enforce the harmonic constraints C = on the boundary and 
then show that this (indirectly) leads to the satisfaction of the Hamiltonian and momentum constraints. 
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An additional complication is that the domain of dependence of the boundary T>(T) by itself is empty. Crudely 
speaking, half of the past directed characteristics from each interior point do not meet the boundary. The Cauchy 
problem is intrinsically coupled with the boundary problem. At the intersection of So and T the Cauchy data 
and boundary data must be prescribed in a consistent way, or otherwise an artificial shock wave will be generated. 
In practice, this is hard to implement without, for example, using an exact solution in the neighborhood of the 
intersection. Moreover, the boundary in general moves relative to the Cauchy hypersurfaces, i.e. the normal to the 
Cauchy hypersurfaces is not in general tangent to T . This complicates the geometric relation between the separate 
3+1 decompositions associated with the boundary and the Cauchy foliation. 

Furthermore, the boundary does not pick out a unique outgoing null direction at a given point (but, instead, 
essentially a half null cone). This complicates the geometric formulation and interpretation of Sommerfeld boundary 
data. This is in addition to the issue of a physical interpretation of the boundary data. Unless the boundary is 
defined by a compact matter distribution, its very existence is physically artificial. This is the situation in numerical 
relativity where a finite outer boundary is typically introduced even though the most important numerical output 
might be the extraction of the gravitational waves that propagate to infinity. The treatment of such an artificial 
boundary can introduce spurious physical effects on the extracted waveform, similar to the effects arising from initial 
Cauchy data which contains spurious gravitational waves. Here I concentrate on the geometrical aspects of the 
boundary treatment but the underlying methods can also be used to improve the physical properties of the treatment 
of an isolated gravitating system e.g. by the construction of boundary conditions which lead to asymptotically 
vanishing reflection coefficients from a sufficiently round boundary for increasingly large radius. 



III. STRONGLY WELL-POSED CONSTRAINT-PRESERVING IBVP WITH SOMMERFELD 

BOUNDARY CONDITIONS 



I begin the discussion of the geometrization of the boundary data with the analytic formulation of a strongly well- 
posed treatment of the IBVP for the Einstein equations in the harmonic gauge based upon 0, [l2| . The harmonic 
coordinates x a — (t,x l ) induce a foliation of the boundary by t — const surfaces Bt- This foliation determines a 
unique future-timelike vector unit T p which is orthogonal to B t and tangent to the boundary. Here the metric has the 
status of an unknown to be determined by solving (| 1 . 5[) . so that other metric related quantities, such as the norm of 
T M , have similar status. Along with the unit spacclikc normal N p which points outward from the boundary, this leads 
to an orthonormal tetrad (T M , N^, Q M .Q M ) at each point on the boundary, where Q p is complex null vector tangent 
to Bt with normalization 

Q M Q M = 2, Q^Q M = 0. (3.1) 

The tetrad is uniquely determined up to the spin-rotation freedom Q p — > e %s Q^ . Uniquely associated with this tetrad 
(independent of the choice of Q M ) are the outgoing and ingoing null vector fields if = T p + N p and L p = — N p , 
respectively. 

The outgoing null vector K p allows us to pose Sommerfeld boundary conditions. In :3j it was shown that a properly 
constructed hierarchy of Sommerfeld boundary conditions of the type 

K^df,Y a = Sommerfeld data, 7^ = y/^gg^ (3.2) 

leads to a strongly well-posed IBVP. In addition to the continuous dependence of the solution on the initial Cauchy 
data, strong well-posedness implies boundary stability, i.e. that the solution for the metric and its derivatives can be 
estimated in terms of the boundary data [131 ] . 

Certain components of the Sommerfeld data are unconstrained. These unconstrained data are picked out by the 
projection tensor Q 



which projects a 1-form V M into the (-£T M , Q M , Q^) subspace. 

The freely prescribed Sommerfeld data q p(T are the 6 components 

ppp$K»d l tf*f i = q p °, (3.3) 

where q p ° ' L a = 0. The remaining boundary conditions enforce the harmonic constraints C u \t = by expressing them 
in the form 

= -d p Y" = +h p K p d^ - P^d^ ~ P^K^d^ = 0, (3.4) 
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which provide Sommerfeld conditions for the remaining L p components of K p d p j p ^. 

Strong well-posedness results from the hierarchical structure of the boundary conditions (|3.3[) and (|3.4[) , i.e. the 
boundary conditions must form a sequence whose Sommerfeld data only depends on prior members. An example, 
which plays an important role in Sec. IIV1 is to prescribe the unconstrained components (|3.3[) in the sequence 



-K p K a K»d^<>° = -V~gqKK, (3.5) 

{QpKvK* -^KpK^dyrf =- V^gqxQ, (3.6) 

(Q P Q„K» -^K p K a L^d^" a = -V~gq Q Q, (3.7) 

{^QpQcrK^ — Q P K a Q^)d^ prJ = -V^gqQQ (3.8) 
and then to add on the constrained components (|3.4[) in the sequence 

2V=S^r" = {K p L a K p + K p K a L p - K p Q a Q^ - K p Q a Q^) d^ p ° = 0, (3.9) 

2y/=g~Q^ = (L p Q a K^ + K p Q a L p -Q p Q a Q^ + Q p Q a Q p )d^ pa = 0, (3.10) 

2^L p T p = (L p L a K» + K p L a L^ - Q p L a Q p - Q p L a Q p ) d^ pa = 0. (3.11) 

The sequential structure is determined by the order 

(K p K ff ), {KpQo), (QpQa), (QpQa), (K p L ff ), (Q p L a ), (L p L a ) 



in which the components of d p 'y f " T enter into the Sommerfeld boundary condition for K p d p j f " 7 . 

The boundary conditions (|3.5|) - (|3.8|) can be modified by the addition of lower differential order terms without 
affecting strong well-posedness. Such terms can be used to reduce spurious back reflections from the boundary of an 
isolated system. (See [llj for more details.) In addition, the boundary conditions (|3.9p - (|3.11[) can also be modified 
in the form 

-H»(x.g) = 

to include gauge forcing terms H p corresponding to a generalized harmonic formulation. Other examples of Sommer- 
feld boundary conditions that fit into the hierarchical structure and that ensure strong well-posedness are considered 

The physical content of these boundary conditions can be clarified by considering the case of a linearized wave 
incident on a plane boundary. The first three pieces of free boundary data {qxK, qKQi 9qq) are related to the gauge 
freedom, i.e. they can be set to zero without loss of generality. The next piece of data (^qq) controls the gravitational 
radiation. The remaining conditions enforce the constraints C = on the boundary. This is consistent with the 
physical expectation that two pieces (or one complex piece) of data are necessary to describe the two polarization 
degrees of freedom of a gravitational wave. More generally, the boundary has curvature, which necessitates additional 
nonzero boundary data, as discussed out in [Tlj . In general, the curvature of the boundary combines with the radiation 
modes in a way which cannot be cleanly separated. This is one of the complications emphasized in the treatment 
given in 0. It is an issue I return to in Sec. IVIBI I next consider how to reverse engineer the above Sommerfeld 
boundary conditions to clarify their geometrical content. 

IV. BOUNDARY GEOMETRY 

I proceed using 4-dimensional notation to describe the 3-dimensional geometrical objects intrinsic to the boundary. 
At the most primitive level, these are the intrinsic 3-metric 

H pv =g pv ~N p N Vl (4.1) 

its corresponding 3-connection D p and the extrinsic curvature 

= H p V p N„. (4.2) 

However, those objects are insufficient to formulate Sommerfeld boundary conditions because they do not determine 
a unique outward null direction. 
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For that purpose, we introduce an evolution gauge vector field f tangent to the boundary. This is in analogy with 
what is done in a numerical evolution the Cauchy problem but for the moment we need only define i M on the boundary 
and not in the interior. It is not necessary that t p be timelike but its flow must map, in a future-directed sense, the 
edge Bq (where the boundary intersects the initial Cauchy hypersurface) into a smooth foliation Bt of the boundary. 
The foliation is parameterized by a time function t satisfying £ t t = 1, with t = on Bq. (Here Ct represents the Lie 
derivative with respect to In turn, along with the intrinsic metric of the boundary, the 2-surfaces Bt determine a 
unit timelike normal field tangent to the boundary according to — —AD p t = — Aff^V^t. The corresponding lapse 
A and shift intrinsic to the boundary are defined according to = AT^ + . 

The intrinsic 2-metric of the boundary foliation is given by 

Q [iv Q{ P Qv) H plJ T^jTu 

where we have again introduced a complex null vector with the normalization (|3.1| . The extrinsic curvature of 
this foliation associated with the normal T M is 

This structure is sufficient to geometrize the last piece of free Sommerfeld data qQQ in the hierarchy (|3.5|) - (|3.8j) 
by relating it to the extrinsic curvature tensors of B t according to 

Q p Q a (ICp a + K pa ) = ^Q p Q a K»d p g pa - Q p K a Q»d p g pa = q QQ . 
Equivalently, we can express this data in terms of the optical shear a of the outgoing null congruence determined by 

2<j = Q p Q°V p K a = q QQ . (4.3) 

In linearized theory, a is the complex data for an ingoing gravitational wave incident on a plane or spherical boundary. 

The free data for the first 3 Sommerfeld conditions (|3 . 5[) - (|3.7[) cannot be similarly expressed in terms of the 
geometry of the boundary and its tangent space. For example, the first condition, rewritten in terms of the metric, is 

X -K p K a K p d p g pa =q KK . 

Since = T M + N^, this contains pieces such as N p N a N^d p g pa which can be set to zero by a diffcomorphism 
leaving the boundary and its extrinsic curvature intact, e.g. by introducing Gaussian normal coordinates adapted to 
the boundary. A geometric formulation of these boundary conditions requires additional structure to eliminate the 
effect of this gauge freedom. One convenient approach is to introduce a preferred background geometry, as described 
next. 



V. INTRODUCTION OF A BACKGROUND METRIC ASSOCIATED WITH THE INITIAL DATA 

Geometrization of the boundary conditions (|3.5p - (|3.7|) can be accomplished by introducing a background metric 

g pv with its associated Christoffel symbols and curvature tensor R P<J ( P - Setting f pi/ = g pv ~ g^ u , the gauge 
freedom in first derivatives of the metric is then fixed relative to the background by the tensor field 

(In order to avoid confusion, raising and lowering of indices will only be done with the physical metric g pv unless 
otherwise noted.) 

The harmonic constraints are now modified to the covariant form 

c p ■■= <T (r£„ - v%) = o, 

or 

C ■■= g pv (r^ - f £„) - H" = o 
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if a gauge forcing term H p (x,g) is included. 

Einstein's equations for reduce to the quasilinear wave system 

<? P "V P V CT /^ = 2gxrg pa C\ p C\ a + AC P a ^g v) xC x pT g aT - 2g^ R X Mf ,g v)x + 2M (p H iy) (5.1) 

for which the formulation of a well-posed IBVP goes through exactly as before. (See [ll| for details.) The null tetrad 
(K p , L M , Q p , Q p ) associated with the physical metric g^ is again defined on the boundary in terms of the foliation 
constructed by introducing an evolution field tangent to the boundary. This null tetrad is then used to prescribe a 
hierarchical set of Sommerfeld boundary conditions. In terms of the analogue of ()3.5() - (|3.11|) take the hierarchical 
form 







l K p K°K p V p f pa 


= QK, 


(5.2) 




{Q p K a K p \ - 


^K p K°Q p \)f pa 


= Qq, 


(5.3) 




(L p K a K p "\7 p - 


l -K p K°L»V p )f pa 


= 1L, 


(5.4) 




(IqpQ'K^ 


- Q p K°Q p %)f pa 


= iQQi 


(5.5) 


(Q p Q a K p + 


K p K a L p - K p Q a Q fl - 


- K p Q°Q p ) \f pa 


= -2K p H pi 


(5.6) 


(L p Q a K p - 


YK p Q a L p - K p L a Q p 


+ Q p Q a Q p ) V p f pa 


= -2Q»H^ 


(5.7) 


{L p L a K p 


+ Q p Q a U l - Q p L a Q» 


- Q p L°Q p ) %f pr7 


= -2L P H^. 


(5.8) 



Note that these conditions can be modified by taking linear combinations consistent with the hierarchy. For example, 
by subtracting (|5.2p we can replace (|5.4[) by 

(N p K a K^ - X -K p K a N^)] pa = q N . (5.9) 

The Sommerfeld boundary conditions (|5.2|) - ()5.8|) for j pa are now in covariant form and the geometrization the 
boundary conditions reduces to eliminating the ambiguity in the choice of the background metric. We do this by 
extending the evolution vector field to a neighborhood of the boundary. This fixes an "evolution gauge" which 
allows us to Lie transport the initial Cauchy data into the future. 

The simplest choice is of background metric is to set 

V, I t=Q = fl7w|t=0 j ^tP^ = 0. (5.10) 
However, higher order compatibility with the initial Cauchy data can be incorporated by setting 

9^=0=9^=0, ^tQfiu = £tSVlt=o- (5-11) 

A further option of (|5.11l) is to to require that the the background metric is harmonic. This is most readily expressed 
in the coordinates x p — (t, x l ) adopted to the evolution according to 

Ctt = l, £ t ^ = 0. (5.12) 

Then the densitized background metric 

J^~gf v = + b pv t + -c^t 2 + -d pu t 3 , (5.13) 
2 6 

where a p,v and b lJ are determined by the initial Cauchy data, satisfies the harmonicity condition provided 

c w = , d pi = , 

b u = -d,a u , b tl = -dja ij , 

c" = -dib u = aidjoV , 

c li = -djb lj d u = -d lC tl = ;),<),!>■ 1 . 

(Here g^ , the inverse of is used to raise indices.) Other background metrics associated with Lie transport of the 
initial Cauchy data could also be used. 
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For both the alternatives (|5. 1 If) and ( I5.13[) , = <7 M „ — </„„ has homogeneous (vanishing) initial Cauchy data 
so that they provide the first approximation to a standard iterative scheme to prove existence of solutions of the 
quasilinear problem (l3| . From a practical point of view, they also automatically ensure a C 1 compatibility between 
the initial data and boundary data, and thus in a numerical evolution would reduce spurious high frequency waves 
arsing at the edge Bo- Such a homogeneous version of the initial value problem can also be useful in numerical 
simulations by eliminating the effect of extreme nonlinear behavior. 

In all three cases, (|5.10[) . (15. lip and (|5. 13[) . the background g is uniquely defined by the choice of the evolution 
vector t M . This applies to any metric based formulation of the Cauchy problem for Einstein's equations. Note that 
all such choices of i M are related by diffeomorphism, so although tP is a critical ingredient it is a pure gauge object. 
It does not determine metric information without further knowledge of the lapse and shift. 

Given a choice of evolution vector tP, the diffeomorphism freedom in the solution to the harmonic IB VP can 
be completely eliminated by tying it to the adapted coordinates determined by a specific coordinatization of So 
according to (|5.12[) . Then the results of ll| guarantee the local existence of a unique solution to the generalized 
harmonic equations ()5.1|) subject to the boundary conditions (|5.2|) - ((57 



VI. GEOMETRICAL INTERPRETATION OF THE SOMMERFELD DATA 

A geometrical interpretation of the Sommerfeld boundary data in (|5.2p - (|5.5p can be made by replacing the 

derivatives V M / pcr in terms of derivatives of R a . The identity (V p - V p )R a = -C^K^ = \RP - V p f pa - V a f pp ) 
can be used to rewrite (|5.5p in terms of the optical shear of the boundary foliation, 

2{a- a) = Q p Q a (V p - V p )R a = \f> p Q° K^ p j pa - Q p R a QP\f pa . = g QQ . (6.1) 

Analogous to (|4.3p . this equates the Sommerfeld data <?qq to the shear of the outgoing null hypersurfaces determined 
by the boundary foliation, relative to the background value of the shear. Alternatively, the spin-weight dependence 
of this relative shear can be removed by expressing it as the rank two symmetric tensor 

fa ■= 2(a^ - a pv ) = (QP^QZ - iQ^XVp - V p )K a . (6.2) 

The identity 

V p K p (V p - V p )RP = (V P R°RP - ^K p K°VP)V p f pa , 

which holds for any vector field V M , can be used to interpret the boundary data for the leading conditions (|5.2p - (|5.4p 
in the hierarchy. The geodesic acceleration of RP with respect to the physical metric and the background metric are 
a' 1 = R"\J V RP and d M = R V \J V RP. Thus (57J) - (53) can be re-expressed as 

K p (a p - &>) = ^K p K°K p V p f pa = q K , (6.3) 
Q p (aP - d") = (Q p K a K p - ^KPK°QP)V p f pa = q Q , (6.4) 
L p (aP - dP) = (LPK°KP - ^KPK°LP)V p f pa = q L , (6.5) 

where the Sommerfeld data {qx, qQ, Ql) is related to the acceleration by the components of 

qP = K " <*")■ (6.6) 

Note that neither a p nor are uniquely determined by our prior construction which only defines RP in the tangent 
space of the boundary. However, their difference is geometrically and uniquely defined, even though either aP or af 1 
can be individually set to zero by extending RP inside the boundary as the tangent vector to the affinely parametrized 
null geodesic determined by either the connection V M or respectively. 

Equations (|6.ip - (|6.5p show how the unconstrained Sommerfeld data can be described in terms of geometrical objects 
consisting of the shear and acceleration of the outgoing null vector RP relative to their background values. This raises 
the question whether another geometrical description can be given which does not introduce a background metric. For 
the shear a, we already saw that this was possible by using (|4.3p . We can investigate the case of a M by decomposing 
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it into the pieces a M = T U V ' V K^ + N VJ \7 V K^. The first piece T V V U K^ is intrinsic to the boundary geometry and can 
be expressed in terms of the extrinsic curvature of the boundary foliation according to the components 

K lt T v V v Ki* = 0, 

Q^VvK" = Q^TK^ + Q^logA, (6.7) 

where A is the lapse intrinsic to the boundary and A^ t „ is its extrinsic curvature. Consequently, if we were to 
prescribe Dirichlet boundary conditions for the metric then the terms N^V^K^ would not enter and the data would 
be "boundary intrinsic" . No background metric would then be necessary. However, the strong well-posedness of the 
IBVP has not been established for the case of Dirichlet boundary conditions. Furthermore, from a practical viewpoint, 
homogeneous Dirichlet boundary conditions are of the reflecting type so that, with the correct prescription of Dirichlet 
boundary data, a gravitational wave could propagate across the boundary but in numerical applications the error 
would be trapped in the grid. 

It is the N V V piece of a M which requires an extension of K M inside the boundary. This leads to a more compli- 
cated geometric description, which depends upon the particular way that is extended. For instance, a boundary 
defining function $ might be introduced such that $|r = and N^\r = r]V^\r, with rj — (g^ V^V^) -1 / 2 . This 
defines an extension of = ryV M $ to a neighborhood of the boundary, where rj plays the role of the "lapse" for the 
boundary defining function. It then follows that 

N^VvN^ = N„N U V V log 77 - V M log r? = log r,. (6.8) 

Next, by extending off the boundary by the parallel transport N V V = 0, we have N v V y K^ = N v V y N^. 
Thus, by putting together the pieces (|6 . T|) and (|6.8[) . a M can be expressed in terms of the extrinsic curvature and 
the lapses A and rj. However, this requires introducing auxiliary quantities, i.e. $ and the extension of the tetrad 
vectors and off the boundary. Consequently, it does not provide any simpler a geometrical description of the 
boundary data than the use of a background metric, except for frame-based formulations, such as [6j, in which the 
tetrad vectors are evolution variables. 



A. Edge data 

The uniqueness of the solution to the harmonic IBVP requires the specification of the initial values of the lapse 
and shift, i.e. the relation rj M = an M + /3 M between the evolution vector £ M and the unit timelike normal to the 
initial Cauchy hypersurface Sq. For the pure Cauchy problem this data is pure gauge information but for the IBVP 
it contains geometric information at the edge where So intersects the boundary T . This edge data is the hyperbolic 
angle sinh O = n^N^ formed by the intersection, which describes the intrinsic motion of the boundary relative to the 
Cauchy foliation. The tangency of the evolution vector i M to the boundary implies that sinh 9 is related to the shift 
by 

smb. 9 = -tf'Wu. (6.9) 
a 

B. Dynamics of the boundary 

The data for a linearized gravitational wave incident on a plane or spherical boundary in a background Minkowski 
space can be prescribed in terms of the shear of the outgoing null hypersurfaces. However, if the boundary has a 
more dynamical behavior, so that its intrinsic metric and extrinsic curvature change dynamically in time, then this 
description breaks down. In the generic case, such a boundary gives rise to a dynamically changing shear even in the 
absence of linearized curvature. Thus the shear must be coupled with data determining the dynamics of the boundary 
in order to unambiguously describe the full physical or geometrical content of the boundary data. 

In the nonlinear theory, the boundary is dynamically traced out by the integral curves of the unit vector . (It is 
also traced out by the integral curves of i M but this description is devoid of metrical content.) Given the initial value 
of T M , which is supplied by the initial Cauchy and edge data, these integral curves can be constructed in principle 
from the geodesic curvature 



A^ = T V V V T^. 
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(The qualification "in principle" is a reminder that the required spacetime metric and connection are unknowns until 
the IBVP is solved.) However, only the normal component of N^A^ enters into determining the dynamics of the 
boundary. We immediately have 

T^A P = 

and, since Ij, = —AD^t (where A is the lapse internal to the boundary and D M its internal connection), we also have 

QuA** = -Q^DuiAD^t) = -Q^AD^Dut = Q^logA 

Thus the components of A^ tangential to the boundary describe the freedom corresponding to the lapse in the foliation 
of the boundary. 

The essential piece of data that controls N^A^ in the boundary conditions (|6.ip - (|6.5|) is L^q p = L fJi (a fJ - — d p ). A 
straightforward calculation gives 

L^q p = -2N^A P + L^N»{V„ - V V )K» - L»T V V U K^ (6.10) 

which is independent of the extension of the tetrad vectors off the boundary. Since the boundary conditions determine 
a unique solution of the IBVP, (|6.10| uniquely determines the boundary values of N^A^. Hence, the dynamical 
properties of the boundary are controlled in the statement of the IBVP by the data L^q^, although in a very implicit 
manner. 

VII. CONSTRAINT PRESERVATION 

The preservation of the harmonic constraints, 

C p :=<T (T^-fP v )-H"(x,g) = 

follows from applying the Bianchi identities to the reduced evolution system. Consider any formulation of Einstein's 
equations for which the reduced evolution equations — take the form 

W v ._ qhv _ v (m C ^) + I g ^ VpC P + A ^ C <y = 0j (7.1) 

for some smooth coefficients A pu (x,g,dg). (This includes constraint modified versions of the generalized harmonic 
system.) The Bianchi identity V^G^" = then implies a homogeneous wave equation for C M , 

V P V P C P + R P C P -2V P {A PP C' J ) = 0. (7.2) 

If the boundary conditions enforce C p \r — and the initial data enforces C p |s = dtC p \s = then the unique solution 
of (|7.2[) is C p = 0. As a result, the Sommerfeld boundary conditions in the geometrical form (|6.2p - (|6.5[) along with 
(15. 6|) - (|5.8p . which enforce C p \r — 0, lead to a well-posed harmonic IBVP in which the constraints C p = are satisfied 
everywhere. In turn, (|7.1j) implies that the Hamiltonian and momentum constraints G^n^ — are also satisfied. 

While these Sommerfeld boundary conditions can be formally applied to any metric formulation of the reduced 
Einstein equations, an independent check is necessary to determine whether the Hamiltonian and momentum con- 
straints are preserved for formulations which do not explicitly contain an evolution system of the form (|7.1[) . An 
important example is the BSSN system [8|, |9( which is widely used in numerical simulations of binary black holes. In 
current numerical work, the boundary conditions for BSSN evolution systems are applied in a naive, homogeneous 
Sommerfeld form to each evolution variable (cf. [Hj])- Constraint preservation does not hold for harmonic evolution 
with these naive boundary conditions and cannot be expected to hold for other systems. The geometric nature of the 
boundary conditions (|6.2|) - (|6.5p suggest that they could be applicable to the BSSN system, although the boundary 
conditions (|5.6p - (|5.8p enforcing the harmonic constraints would undoubtedly need modification. This issue deserves 
further investigation. 

VIII. SUMMARY 

Beginning with an analytic description of a strongly well-posed version of the (generalized) harmonic IBVP, we 
have shown how the boundary conditions on the metric may be expressed in a geometric form. The end result can 
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be summarized as follows. On a manifold M. with boundaries So and T meeting in an edge Bo, we introduce an 
evolution vector field fr 1 whose streamlines are tangent to the boundary and provide a smooth foliation S t of M. 
which intersects the boundary in in a smooth foliation Bt- Although contains no metric information, it supplies 
the essential gauge information to (i) introduce a background metric on M. by the Lie transport of the Cauchy data 
9ij,v\t=o and Ct9fiu\t=Q on So and (ii) introduce a null tetrad on T which is adapted to the foliation Bt- Sommerfeld 
boundary data is then prescribed for the relative acceleration q^, as given in (16. 6|) . and the relative shear q^v, as 
given in (|6.2p . of the outgoing null vector to Bt (relative to their background values). Along with the hyperbolic angle 
characterizing the initial velocity of T relative to So, this data uniquely determines a harmonic spacetime metric 
(locally in time) by solving (|5.ip . 

These boundary conditions have a hierarchical Sommerfeld form which is beneficial for numerical application. Al- 
though formally they can be applied to any metric-based evolution system, an unresolved issue is whether these 
boundary conditions, or some modification, are constraint preserving for hyperbolic reductions of the Einstein equa- 
tions other than harmonic. 

The geometrization of the boundary conditions for the gravitational IBVP has been presented in 4-dimensional 
form in which the data is prescribed in terms of the spacetime metric and its derivatives on iSo and T. A question 
of further geometric interest is whether the data can also be presented in a disembodied 3-dimensional form, as is 
possible for the Cauchy problem by prescribing the initial data h^ v and k^ v intrinsic to the 3-manifold So- Such 
Cauchy data then determines a spacetime metric solving Einstein's equations which is unique up to a diffcomorphism. 
If there is such a disembodied version of the boundary data for the IBVP in terms of the 3-manifolds So and T , and 
their intersection Bo, it does not appear that it can be as simple as in the pure Cauchy problem. In the treatment 
given here, the evolution vector plays key roles in both dealing with the diffeomorphism freedom and describing the 
boundary data in geometric form. The introduction of allows construction of a 4-dimensional background metric 
based upon the Cauchy data. A disembodied version would at the least require introducing on T , along with 
the construction of a 3-dimensional background metric on T . I leave it as an open question whether an equivalent 
3-dimensional version can be given. 
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